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Abstract: In this paper, we consider a de Sitter thick brane model in a pure geometric 
Weyl integrable five-dimensional space-time, which is a generalization of Riemann geome- 
try and is invariant under a so-called Weyl rescaling. We find a solution of this model via 
performing a conformal transformation to map the Weylian structure into a familiar Rie- 
mannian one with a conformal metric. The metric perturbations of the model are discussed. 
For gravitational perturbation, we get the effective modified Poschl- Teller potential in cor- 
responding Schrodinger equation for Kaluza-Klein (KK) modes of the graviton. There is 
only one bound state, which is a normalizable massless zero mode and represents a stable 
4-dimensional graviton. Furthermore, there exists a mass gap between the massless mode 
and continuous KK modes. We also find that the model is stable under the scalar perturba- 
tion in the metric. The correction to the Newtonian potential on the brane is proportional 
to g- 3r / 3 / 2 / / r 2 ) where j3 is the de Sitter parameter of the brane. This is very different from 
the correction caused by a volcano-like effective potential. 
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1. Introduction 

The possibility that our world is trapped in a four-dimensional submanifold (called brane) 
embedded in a fundamental multi-dimensional space-time (called bulk) has been increasing 
interest during recent years (see [|], Q for summary of introduction, and reference therein). 
The strong motivations originate from the string/M theory. With the help of the brane 
scenario, one could possibly solve some disturbing problem of high-energy physics, such 
as the hierarchy problems (the problem of why the electroweak scale Mew ~ 1 TeV is so 
different from the Planck scale M p i ~ 10 16 TeV ) and the cosmological constant problem 
0, [|, [|, ||, [7], |8|. In this scenario, Standard Model matter fields are confined to a 3- 
brane, while the gravity could propagate in the whole bulk. The most well-known models 
within are Arkani-Hamed-Dimopoulos-Dvali (ADD) ||, Q and Randall-Sundrum (RS) |5[ |6[ 
models put forward in the end of 90's. They can achieve a Newtonian gravitation law in 
macroscopic scale, which is compatible with the observational fact. What preeminence in 
RS model is that it suggests a possibility that the extra dimensions may not be compact 
any more, namely the size of extra dimensions can be even infinite. 

However, both the models above are assumed that the brane is infinitely thin. Although 
with the idealized models many interesting results have been obtained, from a more realistic 
point of view, a brane should have the thickness. For this reason, the emphasis of study has 
shifted toward the thick brane scenario |, 0, 0, [12], [D| [14], f|, g§ [17|, [D| [l|, J3J [22 



H H, H, ||, |3, H, |H, H H, H, H . An interesting feature of thick brane scenario based 
on gravity coupled to scalar is that one can achieve a brane naturally without introducing 
a delta function by hand in the action to create a brane. And the scalar could further 
provide the "material" to make the thick brane. In many multidimensional field theories 
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coupled to gravity there are solutions of topological defects. They have lead to a richer 



variety of thick brane worlds [16|. 

Most investigations of brane worlds are considered in Riemann geometry. In this 
paper, we are interested in the thick branes based on a Weyl-integrable geometry, which is 
a generalization of Riemann geometry. 



f 7yo,Q=o 
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As is shown in ([O]) [34], we can see there are two type generalizations of Riemann 
geometry, namely, the Riemann-Cartan and the Weyl geometry. In Riemann-Cartan ge- 
ometry, the connection is asymmetric and the antisymmetric piece of the connection is 
represented as torsion T introduced by Cartan. In contrast to Riemann-Cartan geometry, 
the connection is no longer metric compatible in Weyl geometry, i.e., it involves a geomet- 
ric scalar u) in the definition of connection and the covariant derivative of metric tensor 
is non-vanishing. The field related to this violation of the metricity condition is called 
nonmetricity Q. So Weyl geometry allows for possible variations in the length of vectors 
during parallel transport. When T = the Riemann-Cartan geometry degenerates into 
Riemann, and when Q = the Weyl geometry also degenerates into Riemann geometry. 
More precisely, a Weyl geometry is an affine manifold specified by (<7Afjv 3 w), where gMN 
is the metric tensor and uj is the geometric scalar uj. Since a conformal transformation 
can map a Weyl manifold into a Riemann one, the particular type of gauge geometries is 
called conformally Weyl or Riemann integrable space-time J35|] , On the other hand, Weyl 
integrable manifold is invariant under a Weyl rescaling (see, ( |2.2| ) for detail), when this 
invariance is broken, the Weyl scalar function uj transforms into an observable field which 
generates the smooth thick brane configurations. Thus in Weyl geometry the fundamental 
role in the generation of thick brane configurations is ascribed to the geometric Weyl scalar 
field, which is not a bulk field now. 



In Refs. [14, 35, |36|, 37], the authors have studied brane world scenario in the frame of 
Weyl geometry, and achieved various solutions respected to Minkowski thick branes. For 



most of these branes |35|, ^6[, there exist a single bound state which represents a stable 
4-dimensional graviton and the continuum spectrum of massive modes of KK excitations 



without mass gap for a volcano potential. While in [14], the KK spectrum is quantized for 
an infinite square-like potential wall, and in p7p , there exist one massless bound state, one 
massive KK bound state and the delocalized continuum spectrum for a modified Poschl- 
Teller potential. These give a claim that Weylian structures mimic, classically, quantum 
behavior does not constitute a generic feature of these geometric manifolds pCl]. In Refs. 



1 38, 39], the authors have considered the localization and mass spectrum problems of matter 
fields on these various Minkowski thick branes. An interesting study in |3^] shows that 
for scalars there are two bound states (one is normalizable massless mode), for spin one 
vectors there is only one normalizable massless bound state, and for spin half fermions, the 
total number of bound states is determined by the coupling constant rj: when r/ = 0, there 
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are no any localized fermion KK modes including zero modes for both left and right chiral 
fermions, and when 77 > (77 < 0), the number of bound states of right chiral fermions 
is one less (more) than that of left chiral fermions and in both cases (rj > and rj < 0), 
only one of the zero modes for left chiral fermions and right chiral fermions is bound and 
normalizable. 

We will consider a de Sitter thick brane based on gravity coupled to scalars on a Weyl 
integrable manifold in this paper, and try to find a solution of de Sitter thick brane from 
a pure geometrical Weyl action in five dimensions. The organization of the paper is as 
follows: In Sec. |], we briefly review the Weyl geometry. In Sec. ||, we consider a de Sitter 
thick brane solution on the weyl manifold and find a solution of the model. In Sec. we 
briefly consider the gravitational and scalar perturbations of the structure. In Sec. [|, we 
derive the effective Newtonian potential on the brane. And finally, the conclusion is given 
in Sec. ||. 



2. Weyl geometry 

A non-Riemann generalization of Kaluza-Klein theory is based on the following five-dimensional 

d 5 a 



action [14]: 



s w = f d _^M e -h \ R - 3 £ (Vu,) 2 - W (u,) 
J M w 167rG 5 L 



2.1 



where is a five-dimensional Weyl- integrable manifold specified by the pair (gMN,u), 
9mn is the metric (M,N = 0,1,2,3,5) and u is a scalar function. In this manifold the 
Weylian Ricci tensor reads R MN = T MNA - T^ M N + ^mn T pq ~ T mq F np^ where the 
affine connection of Mf is T MN = { MN } - |(w m5 n + u >N 5 M - gMN^ ,P ) with { MN } the 
Christoffel symbol. The parameter £ is a coupling constant, and U(u) is a self-interaction 
potential for the scalar field. Since the scalar field oo enters in the definition of the affine 
connection of the Weyl manifold, the Weyl action is actually pure geometrical. 
Weyl integrable manifold is invariant under the Weyl rescaling, 

g M N -> V 2 g MN , oj^u + lnn 2 , £ + 8^ In ft 2 ) 2 , (2.2) 

where Q 2 is a smooth function on . In general, this invariance is broken by the self- 
interaction potential U(uj). But from the relation (|2.2|) , we know that U —¥ ft~ 2 ?7 is the 
only transformation preserves such an invariance. Thus the potential will be U(to) = Ae _aJ , 
where A is a coupling constant. When the Weyl invariance is broken, the scalar field could 
transform from a geometric object into an observable degree of freedom which generates 
the smooth thick brane configurations. 

In order to find solutions of the above theory, we shall use the conformal technique to 
map the Weyl frame into the Riemann one. Under the conformal transformation gMN = 
e^SMN, the Weylian affine connection becomes the Christoffel symbol T MN — > { MN }, and 
the Weylian Ricci tensor becomes the Riemannian Ricci tensor. In consequence, via the 
conformal transformation, one recovers the Riemannian structure on the manifold MF'. 
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Now the Weylian action (2.1) is mapped into a Riemann one: 

d 5 xW\§\ 



qR 
^5 



R - 3f (Vw) 2 - 6U(u) 



(2.3) 



where U{oS) = e w ?7(w) and all the hatted magnitudes and operators refer to the Riemann 
frame. Thus, one obtains a theory which describes five-dimensional gravity coupled to a 
scalar field with a self-interaction potential in the Riemann frame. 

By considering the variation of the metric g MN and the scalar function u, we have the 
set of field equations from (12. 



Gmn 



£ duj 



(2.4a) 
(2.4b) 



3. The model 



In Refs. [14, 36, 37], the authors studied thick brane solutions to the theory (p4]), and 
they considered the metric that respects four-dimensional Poincare invariance: 



dsl = e 2A{ - v) r]^dx^dx u + dy 2 , 



(3.1) 



i.e., the Minkowski thick brane solutions in Weyl Geometry, where e 2A ^ is the warp factor 
and y stands for the extra coordinate. The localization and mass spectrum problems of 



matter fields on the Minkowski thick branes were discussed in Refs. [38 



The authors of [jl4j have considered the problem of ^-symmetric manifolds. They 
have chosen £ = —(1 + k)/(4k) and left k as an arbitrary parameter except for k = —4/3. 
The following solution was found: 



M(v) 



2k 



3k 



cosh 



In 



cosh 



3A; 



l-k 



2Xy 



3k 



l-k 



2Xy 



2 

4-3fc 



The authors of [35| have considered another simplified case when k 
as an arbitrary parameter except for £ = —1/16. Their solution is 



1 + 16$ 



In 



y/-8X(l + jjgj 

Cl 



cosh [ci (y - c 2 )] 



V-8A(1 + 16g) 
ci 



2(1+16€) 



cosh [ci (y - c 2 )] 



(3.2a) 

(3.2b) 
-4/3 and left £ 

(3.3a) 
(3.3b) 



Obviously, the first solution (|3.2| ) is symmetric with respect to y coordinate, and the 
second solution is also a symmetric one if we perform a translation of y axis. 
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In this paper, we consider a de Sitter thick brane solution in a pure geometric five- 
dimensional Weyl space-time. The line element for the space-time with planar-parallel 
symmetry is assumed as 

ds 2 = guN<lx M dx N 

= e 2A{v) qpLV dx^dx v + dy 2 

= e 2A ^ (-dt 2 + e^'oVoV) + dy 2 , (3.4) 



where /3 is the de Sitter parameter and related to the 4-dimensional cosmological constant 
of the brane. Under the conformal transformation gMN = e w 9MNi the de Sitter metric 
Q3.4j ) is mapped into 

ds\ = gMNdx M dx N 



c 2 



<v) (- dt 2 + e m dx l dx^ + e-^dy 2 , (3.5) 



where 2a = 2 A — u. 

Now the expressions for Ricci tensor and the scalar curvature in the Riemann frame 

read 

R Q0 = (</' + 4a' 2 + \a'J) e 2A - 3/3 2 , 

4- = _ (</' + 4a' 2 + \a'J) e 2 ( A+ ^ Vij + 3/3 2 e 2 ^, 

£ 55 = _ 4 (a" + a' 2 + iaV), 

R = -4 (2a" + 5a' 2 + a' J) + 12^ 2 e - 2cT , 

where the prime denotes derivative with respect to the fifth coordinate y. The five- 
dimensional stress-energy tensor is given by 

Tmn = Rmn — ~9mnR, (3.7) 

thus, its four-dimensional and five-dimensional components are given through the following 
expressions: 



' oo 



-3(a'' + 2a' 2 + -a'c/)e 2A + 3/3 2 , 



Tij = 3 (a" + 2a' 2 + ±oV) e 2 ^^ - 3/3 2 e 2 ^ % , 



(3.8) 



f 55 = 6a' 2 - 6/3 2 e" 2A . 



For simplicity, using the pair of variables X = uj' and Y = 2 A' , we obtain a set of 
equations from (2^), ( |3.5| ) and (p^): 

X ' + 2XY - -X 2 = I^M e - (3.9a) 
2 £ aw 

y' + 2Y 2 - -XY = | - 417 | + 6/3 2 e- 2A , (3.9b) 
2 \£ au; / 

y 2 - 2XY + (1 - f ) X 2 = -2e" w C/(w) + 4/3 2 e" 2A . (3.9c) 
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Though we have got the field equations by performing a conformal transformation 
9AB = e~ u 'dAB to map the Weyl manifold structure into the Riemann one, the solutions 
corresponding to the Eqs. ( |3.9| ) are harder to get than the similar equations in Ref. [14] 
for the flat brane (just set f3 = in Eqs. (|3.9[)). So, in order to simplify the model, we 
make a coordinate transformation of the form dy = e A ^dz. Then the metric (3.5) will be 
rewrote as: 



ds 



P.2 



,2a(z) 



'q ilv dx il dx v + dz 2 ) = e 2a{ - z) (-dt 2 + e m dx i dx i + dz 



(3.10) 



Now repeating our calculation above, the expressions for the Ricci tensor and the scalar 
curvature corresponding to the metric Q3.10 ) in the Riemann frame read: 



R 00 = a " + 3a' 2 - 3p 2 , 



Rij - 
R = 



-(a" + 3a' 2 -3f3 2 )e 2 ?% j , 
-4a", 

(8a" + Via' 2 - 12f3 2 ) e~ 2CT . 



(3.11) 



The stress-energy tensor is given through the following expressions: 



f 00 = -3{a" + a' 2 -f3 2 ), 
f ij = 3{a" + a' 2 -p 2 ) e 2 ^ VlJ , 
T55 = 6 {a' 2 - p 2 ) . 



(3.12) 



So the set of field equations with the new conformally metric ( 3. 10| ) now read 



fa 

1 dU(uj] 
^ duo 



r'2 



J'2 



a ' 
1 



J 1 R 2 

a - p , 



-2a 



3a' 2 - a" + 3/3 2 ) , 



-2(7 



3a' J + J' 



(3.13a) 
(3.13b) 

(3.13c) 



where the prime denotes derivative with respect to z. If we choose the parameter £ = 1/3, 
and define V(lo) = 3U(uj), we can rewrite Eq. ( |3.13|) as follows: 



12 



dV{u) 
dui 



3(a 



12 



a 



2 



-e 
2 



-2a 



3a' 2 - a" + 3/3 2 ) 



e- 2 ° (3a'u' + 



(3.14a) 
(3.14b) 

(3.14c) 



This set of equations is identical with the field equations of the de Sitter thick brane in the 
Riemann geometry ||, fH 0, [H ||, ||, @5|, @1- Here ' 

we can achieve a ^-symmetric 
solution refer to z coordinate when we consider a symmetric domain wall de Sitter expansion 
in five dimensions for a sine-Gordon potential [47, 48 1 

V(u) = 3U{uj) =pcos 2 (qoj) +v (3.15) 
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with p = 15/3 2 , q = l/\/3, v = — 15/3 2 /2 and /? > 0. The warp factor of the metric ( 3.1C| ) 
on Riemann manifold and the scalar uj are 



e M z ) = sech (2/3 
V3 



w(z) 



arctan [sinh(2/3,z)] . 



(3.16a) 
(3.16b) 
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Figure 1: The shapes of the scalar wmz coordinate. The parameter f3 is set to f3 = 0.9 for the 
dashed line, and (3 = 1 for the thin line. 



U(co) 




U(oj) 







/ ^0 




/ x " " x ?\ 
/ / 4-^ 




l ^ 


\\ 




\ \ 

- \\ 


1 1.5 


/ 1.0 


\\ 
\\ 


// 


V\ 


\ °- 5 


\\ 



-1.0 -0.5 



0.5 1.0 



1.0 -0.5 



0.5 1.0 



CO 



Figure 2: The shapes of the potential U{uj) on the Riemann manifold (left) and the potential U(ui) 
on the Weyl manifold (right). The parameter /3 is set to (3 — 0.9 for the dashed line, and j3 = 1 for 
the thin line. 

From U(lo) = e^U^) and Eq. ( |3.15 ), we have the potential in the original action ( |2.1| ) 
on the Weyl manifold: 



U(oj) = 5(3 2 



cos 



u \ 1 



From the relation e 2cr = e 2 (^~ w ), the warp factor e 2A ^ is given by 



2A(z) _ V3arctan[sinh(2/3z)] 



sech (2(3z) . 



(3.17) 



(3.18) 



The shapes of the scalar uj, the potentials U(u) in the Riemannian structure and U (oj) 
in the Weylian one, the warp factors e 2<J on the Riemann manifold and e 2A on the Weyl 
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z 



Figure 3: The shapes of the warp factor e 2a on Riemann manifold (left) and the warp factor e 2A 
on Weyl manifold (right). The parameter f3 is set to (3 — 0.9 for the dashed line, and /3 = 1 for the 
thin line. 



manifold are shown in Figs. |T], [2] and ^, respectively. We can see that the scalar field 
takes values ±\/37r/4 at z — > ±oo, corresponding to the two minima of the potential with 
cosmological constant A = 0. So the scalar is actually a kink solution, which provides a 
thick brane realization of the brane world as a domain wall in the bulk. In the Riemannian 
structure, the potential U{u) and the warp factor are symmetric, but in the Weylian one, 
they are asymmetric. 

Now from dy = e A ^dz, we can calculate the relation between the original space-time 
coordinate y in metric ( |3.5| ) and z in ( p. 10 ). However, we could only achieve a numerical 
function curve y = y(z) in Fig. |] (left). The scalar field in y coordinate is plotted in Fig. 
| (right). 
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Figure 4: The relation of the coordinates y and z (left) and the scalar field to in y coordinate 
(right). The parameter /? is set to /3 = 0.9 for the dashed line, and j3 = 1 for the thin line. 

We can find that the scalar uj(y) is asymmetric but still a kink- like solution in y 
coordinate. We display the shapes of the warp factors e 2a ^ and e 2A ^ in Fig. || The 
figure shows that the warp factor e 2a ^ is no longer symmetric even in the Riemannian 
structure, though the maximum is still at the origin of the y coordinate. So it will not 
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Figure 5: The shapes of the warp factors e 2<y ^ on the Riemann manifold (left) and e 2A ^ on the 
Weyl manifold (right). The parameter (3 is set to (3 = 0.9 for the dashed line, and (3 = 1 for the 
thin line. 



preserve ^-symmetry in the Riemann manifold in the y coordinate. 

The energy density of the scalar matter on the Riemann manifold is given by the 
null- null component of the stress-energy tensor ( 3.1 2| ) : 



3 (a" + a' 2 



p(z) = Too : 

= 9/3 2 sech 2 (2/3z). 



/3 2 ) 



(3.19) 



The shape of the energy density p is plotted in Fig. [6| in z and y coordinates, respectively. 
The figures clearly show that the energy density p distributes along extra coordinate, so it 
does not stand for a thin brane any more, and 1/(3 plays the role of the brane thickness. The 
energy density has a maximum at the origin of the coordinate, and vanishes asymptotically 
at y = ±oo, so the scalar matter mainly distributes on the brane. The energy density is 
symmetric in z coordinate while asymmetric in y. 
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Figure 6: The shapes of the energy density p in z coordinate (left) and in y coordinate (right) on 
the Riemann manifold. The parameter (3 is set to (3 = 0.9 for the dashed line, and (3 = 1 for the 
thin line. 



4. Metric perturbations 

Firstly, we consider the gravitational perturbation in the metric ( 3.10Q : 
ds\ = e 2a{z) [( gflu + h^ix, z)) dx^dx u + dz 2 ] 

a 2a{z) \_ dt 2 + & m dx i dx i + h r X) z )d X ^dx U + dz 1 



(4.1) 



Here h^ u is the tensor perturbation of the metric, and it satisfies the transverse traceless 
condition g, g{|, 



(4.2) 



where V denotes the covariant derivative with respect to the four-dimensional metric q^ v . 
The equation for h^ v is given by jDI 



0. 



(4.3) 



where □ = q^ v V „Vj,. Here we can define the four-dimensional mass of a KK excitation as 



Uh 



fJtV 



Furthermore, we can decompose L„ in the form 

huv(x,z) = e- 3a/2 e^ u (x)ty(z), 



where satisfies the transverse traceless condition 



e/ = VV = 0. 



Now considering ( |4.4| ) and (4.5), Eq. ( |4.3P can be rewritten as 

+ f^-a" + Jja' 2 ) <Z(z) = m 2 *(z). 
So it turns into a Schrodinger equation form, and the effective potential reads 
V g (z) = \a" + V = h 2 [3 - 7sech 2 ( 2/ 3z)] . 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



The spectrum of the eigenvalue m 2 parameterizes the spectrum of the observed four- 
dimensional graviton masses. And obviously, there exists a zero mode which refers to 
m = 0. From Eq. (4.7) with m = 0, we can easily get 



cq sec. 



h 3 /\2(3z). 



It can easily be seen that this zero mode is normalizable, and cq 



(4.9) 

2/3r(5/4) - g a norma ^_ 



v^r(3/4 



ization constant. So the 4D gravity can be produced by the zero mode 



Tg, y, |j. 
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Figure 7: The shapes of the potential V g (z) (concave) and the KK zero mode ^q(z) (convex). The 
parameter f3 is set to /3 = 0.9 for the dashed line, and /3 = 1 for the thin line. 

The shapes of the potential and the zero mode are plotted in Fig. 0. The figure 
shows that the gravitational zero mode localizes around z = 0, and vanishes at each 
side asymptotically. The effective potential is actually a modified Poschl- Teller potential. 
From Q4.S| ) we known that the potential barriers at each side approach to a maximum 



ymax _ Qj-jZ/^. w hen z — > ±oo, so there will be exist a set of continuous KK modes ^> m (z), 
and there is a mass gap between the massless zero mode and the massive continuous KK 
modes. The existence of the gap is universal in various kinds of de Sitter 3-brane model 
1 13, E^]. It can be shown that there is no any other bound KK mode beside the zero mode 



for the potential (4.8). 



Next, we consider the scalar perturbation. Following the arguments in Refs. [13, 49], 



one can get the scalar perturbations in the metric ( 3.10| ) with the longitudinal gauge: 

ds\ = e 2CT ^ [{g M N + <W) dx M dx N ] 

= e 2a{z) [(1 + 2(f)) dz 2 + (1 + 2V>) q^daPda?] . (4.10) 

The equations of the scalar perturbation in the Schrodinger-like form are read as p9|: 



Su = — (-3^ + 3a» , (4.11a) 

cp + 2ip = 0, (4.11b) 
uj' 

Hx,z) = ^/2$(x,z), (4.11c) 

-$"(x, z) + V s (z)<f>(x, z) = U${x, z), (4.12) 

where □ = g^V^V^, 5oo is the perturbation of background scalar u, and V s {z) is the 
effective potential of the scalar perturbation with the expression 



V s (z) = --a" + + t /^~ + 2[^) (4.13) 



Now we decompose $(x, z) according to ( (4.12|) in the form 

$(x,z)=X(x)f{z). (4.14) 
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Then substitute it into Eq. ( 4.12| ), we get 



BX(x) = m 2 X(x), 
-C{z) + V s {z)f m {z) = m 2 f m (z). 



(4.15a) 
(4.15b) 



where m is the 4-dimensional effective mass of the scalar perturbational field. Now Eq. 
(HUD can be explicitly expressed as 



nx( x ) = <rv»v v x(x) = <r [a^du - T%,dx) x(x) 

= (-d 2 + e-Wfididj - 3/3$) X{x) = m 2 X( 



i.e., 



(4.16) 
(4.17) 
(4.18) 
(4.19a) 

T km {t) + 3f3T km (t) + (e-Wj? + m 2 ) T krn (t) = 0. (4.19b) 
The solution of Eq. ( 4.19a|) is simply achieved as S k (x l ) = ce - lS ij k * xJ ; where c is a constant. 



e 2l3t ($ 2 + 3/3$ + m 2 ) X{x) = 5 ij d l d j X(x). 
So X(x) can be decomposed as 

X(x) = T km (t)S k (x i ). 

Then we get 

ftdidjSkia?) = -k 2 S k (x l ), 



Thus we will focus on Eq. (4.19b). The solution is 



J_(\k\ 



-fit 



(4.20) 



where ci and c 2 are integration constants, ? = y |/3 2 — m 2 , and J is the Bessel function. 
Now we expand T km {t) about the point ^■e~ /3t ~ (t — > oo): 



T fcm (t) = Cl 2-5+* f^e-^ 



3_i 

2 p 



3, i 



+ c 2 2 



-/3t 



2 /3 



5,i' 
2^/3 



\P J \\P J 

c ; e -(§4»-0* + O ( e -(f ^-0*) + ^ e -(I^)* + {e~^ + ^) , (4.21) 



where c' x and c' 2 are constants independent of time. 
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Case 1: 



i( with C > 0. Now (gjjj ) 



When m > |/3 , <r is imaginary. We define ? 
can be rewritten as 

T km (t) = c^d^O* + 4e-(^+«*. (4.22) 

Fig. ||a shows that, in this case, 7fc m (i) vanishes at t — > oo, i.e, the scalar perturbation 
disappears after a period of time. So we can see the structure is stable in this case. 

Case 2: When < m? < jj3 2 , <j is real and satisfy < c < |/3. Eq. ( |4.21[ ) shows 
that Tkm{t) is actually suppressed by i, and vanishes at t — > oo, just as shown in Fig. |8|b. 
So the structure is also stable in this case. 



Re[T km (t)],Im[T km (t)] 
3.0r 



3.0 



2.5 
2.0 
1.5 
1.0 
0.5 
0.0 
-0.5 




io o.o 



Figure 8: The shapes of the function Tk m (t) in case 1 (left) and case 2 (right). In case 1, the 
real part of Tk m (t) is denoted by the dash line, and the imaginary part is by the thin line. The 
parameters are set to d x = 1, c' 2 = 2, and j3 — 1. 




02468 10 2468 10 



Figure 9: The shapes of the function Tj- m (t) in case 3 (left) and case 4 (right). The parameters 
are set to c[ = 1, c' 2 = 2, and (3 = 1. 



Case 3: When m 2 = 0, c = f/3, the high order terms O [e'^'^j and O (V^^) 
in the expanding expression ( 4.21| ) are suppressed by t. While the first term in ( |4.21| ) is 



a 



constant c[ = c\. If c\ = 0, Tfc m (t) also converges to zero when t — > ±oo, so the structure 
is stable. If c\ / 0, 7fc m (i) will converge to the constant c\ (see Fig. ^a), i.e., some slight 
scalar perturbations will always exist but not enormously affect the stability of the brane. 
So in this case we can say that the structure is also stable. 
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Case 4: When m 2 < 0, q is real and satisfy q > |/3, the first term in ( 4.21| ) is 
divergent at t — > oo (see Fig. p|b). This means that in this case the scalar perturbation is 
divergent, and the structure is unstable under the scalar perturbation. 

From above discussion, it is shown that the structure is stable when m 2 > 0, while it 
is unstable when m? < 0. 

We next focus on Eq. ( (4.15b| ), where V s (z) is the effective potential of the scalar 
perturbation. For ( 3.16a| ) and ( |3.16b| ), we have 



V s { Z ) = -"-a" + y* + a<"--"- + 2r-\ -6/3 2 
2 4 oj' u' \ lo 'J 

= [3 + sech 2 (2/3z)] . (4.23) 
From the expression of the above effective scalar perturbational potential, we can see that 

V s (z) 



U) 




Figure 10: The shapes of the potential V s (z). The parameter /? is set to (3 = 0.9 for the dashed 
line, and j3 = 1 for the thin line. 

V s (z) approaches to |/3 2 when z — > ±oo. Furthermore, V s (z) satisfies V s (z) > |/3 2 and is 
always convex and non-negative. Fig. shows the potential with virous cases about the 
parameter (3. From the Schrodinger-like equation (|4.15b ), we can conclude that there are 



only KK modes with eigenvalues m? > |/3 2 . This coincides with case 1. So the structure 
is stable under the scalar perturbations. The explicit argument about the stable problem 
of scalar perturbations could be seen in Ref. [[h]]. 

The coordinate system corresponding to the four-dimensional metric q^ v in ( [3.4] ) is not 
a global one. However, even if we use a global coordinate system, we will find the same 
results as the previous discussion, namely, the scalar perturbation will ultimately vanish, 
approach a constant, and diverge for the cases of m 2 > 0, m 2 = 0, and m 2 < 0, respectively. 
So this consequence results from the property of the four-dimensional background space- 
time. On the other hand, Eq. ( 4.15a| ) is just a Klein-Gordon equation. If the background 



is a flat Minkowski space-time, the solution is just a plane wave which never vanishes 
and diverges. But now the background is a de Sitter space-time. There is a exponential 
expansion factor e 2 ^* before the space components, that means the space will exponentially 
expand with time. As a consequence, the symmetries of a de Sitter space-time are very 
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different from that of a Minkowski one. However, the space is still homogeneous and 
isotropic in de Sitter case, so the solution respect to the Eq. ( 4.19a| ) has the same form as 



in a Minkowski space-time, i.e., Sk(x l ) = ce~ iSi ^ k x3 . But the solution of time component 
will change and we find that the amplitude will ultimately vanish especially for m 2 > in 
de Sitter case. Then from Eqs. ( |422| ) and (f01|), we can see that <3?(x,z) will approach 



to zero. Therefore, from ( 4,11| ), we know that the scalar perturbations </>, t/j and Soj will 
ultimately vanish. 

In what follows, we will address the question of how the perturbations couple to matter 
fields on the brane. However, since the scalar perturbations cj>, ip and 5to ultimately vanish, 
their interaction with matter fields will decouple finally. So we only consider the interaction 
of the tensor perturbations couple to matter fields on the brane, which will lead to the 
effective Newtonian potential on the brane. 

5. The Effective Newtonian Potential 

In Section ||, we have considered the gravitational perturbation and achieved a Schrodinger 
equation ( [4.7[ ). In order to have localized four-dimensional gravity, we should require 
that the KK modes of the Schrodinger equation ( (4~7[ ) don't lead to unacceptably large 
corrections to the Newtonian potential in four-dimensional theory. In the thick brane 
scenario the matter fields in the four-dimensional theory on the brane would be smeared 



over the fifth dimension. For simplicity, like in Refs. [jl|, ||, 50, 51, 52|, we consider the 
gravitational potential between two point-like sources of mass Mi and M<i located at the 
origin of the fifth-dimension, i.e., z = 0. This assumption is justified in case when the 
thickness of the brane is small compared with the bulk curvature. The effective potential 
between these two particles is given by the exchange of the massless zero mode and KK 
massive modes. And the zero mode will cause a four-dimensional Newtonian interaction 
potential, the continuum KK modes will produce the correction to the potential. Thus the 



effective potential is given by [51 



mo 



where Gjy = M^ 2 /IQtt is the four-dimensional coupling constant, i.e., the Newton's grav- 
itational constant, M* = {lQ-nG^)~ 1 ^ is the fundamental five-dimensional Planck scale, 
and mo is the minimal eigenvalue at which the continuum KK modes start. In our model, 
mo is n ot zero but 3/3/2 because there is a gap in the mass spectrum. 

As in ||, consider the four-dimensional perturbational metric factors in ( p. 10 ) 



dsl = e 2 ^ [(</> + h^{x)) dx»dx u + dz 2 ] , (5.2) 



We decompose the five-dimensional action (2J3) into a four dimensional and higher dimen- 
sional parts, i.e., 



D Mf J d 5 x^R D Ml J dze 3a ^ J d 4 x/§WR^ = M 2 J d A x^)R^\ (5.3) 
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where is the four-dimensional Ricci scalar made out of g$ = Q^v + h tiu {x). So we can 
read off the relation of the effective four-dimensional Plank scale M4 and the fundamental 
five-dimensional Planck scale M#, 

Ml = Ml £° dze^ = r 2 (^) ^ « 1.2^. (5.4) 

We note that since M4 = M p ; and the fundamental five-dimensional Planck scale M* is 
not fixed, the order of M* is decided by the de Sitter parameter f3. However, an analogue 
relation can be found in RS2 model % g, M| = Mj*/k, where fc appearing in the warp 
factor of RS metric is a scale of order the Planck scale and represents the anti-de Sitter 
radius. Thus M* ~ M v \ in RS2 model. From Eq. (5.4), if we chose /3 ~ M p i, M* is of 
order the Planck scale in our model. 



In order to get the continuum modes, substituting (|4.8|) into (4.7) and introducing a 



new variable I = 2f3z, we can reform the Schrodinger equation as a simple form: 

_*»m _ ^ sec h 2 (Z)^(0 = M 2 f (/), (5.5) 
16 

where the prime denotes derivative with respect to the coordinate I, and the new eigenvalue 
is M = J jj^z — Since in our model, the continuum KK modes start at m > 3/3/2, so 
when we just consider these continuum modes, M will be a nonnegative real number. 
The solution of this equation is given by a linear combination of the associated Legendre 
functions: 

m M (l) = d P(-,iM,tanh(0) + C 2 Q(-,iM,tanh(/)), (5.6) 

where C\, C 2 are M-dependent parameters. When I — > 00, each KK mode wave function 
will approach to a plane wave: 

*m(0 =Ki e iMl + K 2 e~ iMl (5.7) 

with the parameters K\ and K 2 given by 

2Ci - mC 2 coth(vrM) 
K 1 = - 



K 



2r(l - iM) 

C 2 



2 



2C M [2sinh 2 (M7r) + isinh(2vrM)] ' 



where C M = £372 (3 + 8iM + 16M 2 )r(l + zM)r(-§ - 2iM). Further, we chose the 
parameters K\ = K 2 = 1/2 to normalize the plane wave function ^ m (I), so the parameters 
C\ and C 2 are fixed as: 

d = |C M [2cosh 2 (vrM) - i sinh(2vrM)] + -r(l - iM), 
C 2 = -C M [2sinh 2 (Mvr) + i sinh(2vrM)] . 

Substitute C\ and C 2 into Eq. ([5l]), we will achieve the whole wave functions of 



the continuum KK modes *m(0- Now replace I with 2/3 z and M with */ ™ — ^ in the 
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4^2 16 
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Figure 11: The shape of |^E' m (0)| 2 as a function of m. The black dashed line represents |4' m (0)| 2 , 
the red dashed line represents the constant 0.070, and the blue thin line represents the constant 1. 
The parameter is set to (3 — 1. 



function we transform ^m(1) back into *$> m (z) which refers to the KK modes wave 

function in original Schrodinger equation (4.7). We display the curve of |\& m (0)| 2 as a 



function of m in Fig. 11. It is shown that | x I / m (0)| 2 will approach to 1 with the increase 
of the mass m. The reason is that the wave function ty m (z) approximates to a normalized 



2 is a 



plane wave as the eigenvalue m become large. On the other hand, the value of \^ mo 
constant 7r/[r(l/8)r(ll/8)] 2 = 0.070 and independent of the parameter (3. 

Now we focus on the Newtonian potential correction term in Eq. (|T 
integrand as: 

I(m) = ^H* m (0)| 2 . (5.8) 



We write the 



Because it can not be analytically integrated, we will deal with the complex integrand I{m) 
under some simple cases for which the integral can easily achieved. 

Case 1: We replace all the KK mode wave functions \P m (0) with the lowest mode 
* mo (0). In this case, |^ m (0)| 2 will be replaced by vr/[r(l/8)r(ll/8)] 2 = 0.070. Now the 
integrant is simply given by 

p—mT 

(5.9) 



h{m) = 0.070- 

v 

So in this case the correction of the Newtonian potential is 

e -3r/3/2 MlM2 



A[/i(r) = 0.070- 



Af3 



(5.10) 



Case 2: We replace |^ m (0)| 2 with the constant 1, which refers to the limit | ^00(0)! 2 
when m — > 00. In this case, all KK modes are approximately considered as the plane waves. 
So we have 

e -mr 

h(m) = . 

r 

The correction of the Newtonian potential is given by 

g -3r/3/2 MlM2 



(5.11) 



AU 2 (r) 



MS 



(5.12) 
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Case 3: The wave function | x I' m (0)| 2 approaches to a constant 1 when m is large and 
the exponent e~ mr is strongly suppressed by large m, so the integrant is mostly determined 
by the contribution of smaller mass modes. For this point, we can expand | l I / m (0)| 2 about 
the point mo = 3/3/2: 



m 



m 



|<lr m (0)| 2 = 0.037 - 0.079-2- + °- 042 ^4 + ° 



P 2 



4 



m 



fj 6 



So the integrant now can be expressed as 
J 3 (m) 

In this case, the correction can be write as 



m 



0.037 - 0.079-2- + 0.042' 

fj fj 



m 4 \ e mr 



(5.13) 



(5.14) 



AU 3 (r) 



e -3r/3/2 MiM2 



0.072 0.330 0.976 1.512 

O ~r\ Q "t" A nO ~T~ E nQ ~T~ 



/?r 3 



r 4 /3 2 r 5 /3 3 r 6 /3 



1.008 



(5.15) 



In order to distinguish which case above is more accurate compared to the original 



correction of Newtonian potential, we can compare the precise integrand (5.8) with the 
above three approximate cases (|5.9| ), (|5.11 ) and ( 5.14 ) at some fixed distances r of the two 
particles. Since the area under an integrand curve represents the integrate value of it, so by 
comparing the areas of three cases with the precise condition, we can fixed a suitable form 
of the correction. Now we display the curves of the precise and above three approximate 



cases with some different fixed distances in Fig. 12 and 13. 




10 



Figure 12: The integrand curves of the precise and three approximate cases with fixed distances 
r = 0.8 (left) and r = 1 (right). The black dashed thick lines represent the precise integrand curve, 
the red dashed thin lines represent case 1, the blue thin lines represent case 2, and the orange thick 
lines represent case 3. The parameter is set to /3 = 1. 



We find that the integrand ( |5.14 ) in case 3 largely deviates from the precise case when 



r is small (Fig. 12) but well matches with the precise case when the distance r is large 
enougn (Fig. p|). The reason is that we replace the function |^ m (0)| 2 in the integrand 
J5.8|) with a polynomial of m. The polynomial is a good approximation at small m, and 
even at large m, as long as the distance r is also large. Large r will lead to the exponent 
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I(m) l(m) 




Figure 13: The integrand curves of the precise and three approximate cases with fixed distances 
r = 5 (left) and r — 10 (right). The black dashed thick lines represent the precise integrand curve, 
the red dashed thin lines represent case 1, the blue thin lines represent case 2, and the orange thick 
lines represent case 3. The parameter is set to /3 = 1. 



e -mr g rea ^iy suppressed, so the polynomial is still a good approximation. However, when 
r is small enough, e~ mr will not be very little, but the contribution of the polynomial with 
large m will be great, in consequence there will be a obvious derivation with the precise 
case. This is also can be seen in Eq. ( |5.15| ), when r is small enough, the main terms 
that will contribute to the correction are the high order terms of 1/r 2 , but when r is large 
the main term will be 1/r 2 . So for this consideration we will not use the correction of 
case 3 as the approximation of the final correction. On the other hand, we can read from 
the two figures that the precise integrand curve is always between the curves of case 1 



and case 2, and this is obviously, since the Fig. 11 shows that 0.070 < | x I' m (0)| 2 < 1, so 
these integrands have the relation I±(m) < I(m) < I^im). Hence, the precise correction of 
Newtonian potential AU(r) must satisfy AU\{r) < AU(r) < AU2(r), i.e., 

-3r/3/2 M M e -3r/3/2 M M 

0-070—-^ < A U(r) < V^- (5-16) 

M| r 2 v ; M| r 2 y ' 

Though the above discussion, we obtain the suitable correction of Newtonian potential: 

AZ7(r) = k—, (5.17) 

where the parameter k is a constant satisfied 0.070 < k < 1. 

So in our mode the effective Newtonian potential can be written as 

U{r) = G N —^ + K—z ^- . (5.18) 

This is greatly different from the correction caused by a volcano-like effective potential 



1 51, 52]. We can find when the distance r of the two particles is large, the 1/r 2 term is a 
high order term compared with 1/r, and the exponent e _3r ^ //2 is also greatly suppressed. 
So in this situation, the correction is tininess and can be negligible. However, when r is 
small, e _3r ' 3 / 2 /r 2 will be large and could be the main term of the potential. 
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6. Conclusion 



In this paper, we consider the de Sitter thick brane solution in Weyl geometry. By per- 
forming the conformal transformation to map the Weyl structure into the familiar Riemann 
one, and further, via a coordinate transformation dy = e A ^dz to transform the metric 
into a conformal one, we transform the structure to a de Sitter thick brane on Riemann 
manifold. The solution is 



U(u) = 5/3 2 



cos 



u \ 1 



VsJ 2 
3inh(l 

2A(z) _ e v / 3arctan[sinh(2/32)] sec j 1 (2/3z) . 



s/3 

uj(z) = —— arctan [sinh(2/3 



(6.1) 

(6.2) 
(6.3) 



The scalar field a; is a kink solution and preserves symmetry with respect to z coordinate, 
while the scalar potential U (lj) and warp factor e 2A are asymmetric. But actually, z 
coordinate is not a real physical one. After transforming back to the physical y coordinate 
with dy = e A ^dz, we find that all magnitudes are not preserving symmetry with extra 
coordinate y. The energy density is still localized around the origin of coordinate, while 
our matter distribution on the de Sitter brane is asymmetric along the extra dimension. 

Then we consider the perturbations in the metric ( |3.10| ) in Riemann structure. For 
gravitational perturbation, we get a Poschl- Teller-like potential for the gravitational KK 
modes. It is found that there is a normalizable massless gravitational zero mode localized 
on the brane, and exists a mass gap between the zero mode and the massive continuous 
KK modes. The existence of such a mass gap also conforms to the universal phenomenon 
in various de Sitter 3-brane models. And for scalar perturbation, it is shown that the scalar 
perturbations will ultimately vanish, so the solution is stable in this perturbation. 

Finally, we discuss the effective Newtonian potential on the brane, and find the cor- 
rection term is proportional to e - 3r ^/ 2 / r 2 ) which is greatly different from the correction 
caused by a volcano-like effective potential. This result shows that the the perturbations 
don't lead to unacceptably large corrections to the Newtonian potential in four-dimensional 
theory. 
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